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Abstract 

The horizontal lateral density of a cosmic air shower with a non¬ 
zero zenith angle is asymmetric. The asymmetry consist of a stretching 
of the iso-density contours to ellipses and to a shift of the center of the 
elliptic contours with respect to the core of the shower. The shift is 
caused by atmospheric attenuation. The modeling of the attenuation 
results in an equation for the shift as a function of zenith angle and the 
size of the iso-density contours. A more accurate equation is obtained 
by investigating the shift in lateral densities of simulated showers. It is 
shown how the shift can be incorporated in an elliptic lateral density 
function. A linear approximation for the shift allows for an analytical 
solution for the shifted elliptic density. Its predictions for the polar 
variations of the density are compared with data of simulated showers. 
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1 Introduction 


A lateral density function (LDF) describes the density as a function of the 
radius with respect to the core of a shower. For vertical air showers the 
horizontal plane coincides with the plane of the front of the shower and the 
iso-density contours in the horizontal plane are circles. A polar symmetric 
LDF is of application for vertical showers and for polar averaged densities 
of inclined showers. For inclined showers the iso-density contours are rather 
ellipses mm- As known, the centers of the elliptic iso-density contours do 
not coincide with the shower core, see Fig. 



Figure 1: Impression of the lateral density by means of iso-density contours 
and the different positions of the shower core and the point of intersection 
of the major axis and the minor axis, the center, of an iso-density contour. 


The distance between the shower core and the center of an elliptic contour 
will be denoted as the ‘shift’. The application of an elliptic LDF instead of 
a polar symmetric LDF increases the accuracy of the reconstruction of an 
inclined air shower observed with detectors in a horizontal plane. The accu¬ 
racy of the reconstruction of an inclined air shower can be increased further 
if the shift is taken into account. In order to give already an impression an 
LDF-A solely based on the projection, thus without a shift, and an LDF-A 
including the shift, are both plotted for the polar density for an average 100 
PeV shower with zenith angle 45° at a distance 100 m from the core in Fig. 
1^ We see the additional angular density variations caused by the shift are 
of the same order as the angular density variations caused by the projection. 
This suggest that if the ellipticity is taken into account for reconstruction 
purposes, then the shift might be taken into consideration as well. 
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Figure 2: The polar density according to projection with a shift (solid) and 
projection without a shift (dashed) for an arbitrary 100 PeV shower with 
zenith angle 45° at a distance 100 m from the core. The horizontal line 
(dotted) is the mean density. 


The main purpose of the paper is the construction of an asymmetric density 
function which includes the shift. To this end the shift will be investigated 
for different primary energies and different zenith angles. The shift is caused 
by the attenuation of the shower. For the electron part this is the atmo¬ 
spheric attenuation which can be modeled to a certain extent. For the muon 
part the attenuation is mainly due to decay. The decay of the muons can be 
modeled. However, a decaying muon contributes an electron to the electron 
part. The two distribution are therefore intertwined. It therefore does not 
make much sense to consider the shift for electrons and muons separately. 
Besides it would require to consider the ratio of the densities of the elec¬ 
trons and muons in a model. Although we are interested in the shift of the 
combined density of electron and muons together, we will for the model re¬ 
strict to atmospheric attenuation. Even with this restriction the accuracy 
of the model is limited for several reasons of which the ignorance of local 
attenuation is the most important. Nevertheless, the model result gives an 
indication of the way the shift depends on zenith angle and on distance to the 
core. Accurate values for the shift are determined from lateral densities of 
MC showers. In plots of the determined shifts the model prediction will still 
be plotted for reasons of comparison. Furthermore it will be shown how the 
asymmetric LDF including the shift is constructed from a polar symmetric 
LDF. To avoid length we will denote a polar symmetric LDF just as LDF 
and an asymmetric LDF as LDF-A. For a clear distinction we will denote the 
density in the front plane as p and the asymmetric density in the horizontal 
plane as v. 

The contents of the paper can be divided in three parts: the modeling of 
the shift (Section 2-5), the determination of the shift (Section 6-8) and 
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the construction of a polar density function including the shift (Section 9 - 
10). In Section 2 we consider a cylinder model for the shower in a suitable 
coordinate system. The consequences of the cylindrical projection will be 
considered for the situation without and with a shift. In Section 3 we model 
the effect of atmospheric attenuation of the shower on the lateral density. 
An analytical approximation for the shift is derived in Section 4. In Section 
5 we obtain a comparable result on the basis of a cone model for the shower. 
In Section 6 a description is given of the method used for the investigation of 
the shift on the basis of horizontal densities of simulated showers. In Section 
7 some general results will be presented as obtained from simulated showers. 
In Section 8 we focus on the behavior of the shift for the combined density 
of electrons and muons together. The muon energy deposit in scintillators 
is practically similar to the electron energy deposit |3]. A combined lateral 
density of electrons and muons is therefore of interest for scintillator based 
observatories. The detecting efficiency of scintillator detectors become small 
for densities below 0.5 We therefore will focus on combined densities 

larger than 0.5 m“^. We will see that in this region the shift is independent of 
shower size. Moreover, the relation between shift and the size of the elliptic 
contour is almost linear. The shift of the combined density will be compared 
with a proposed linear approximation. In Section 9 the polar density is con¬ 
sidered. It is shown how to convert an LDF to the LDF-A including the 
shift as a function of radius r and polar angle a and parameterized by the 
zenith angle 9. For the outline of the procedure we conveniently restrict to a 
proposed linear approximation for the shift since it allows for an analytical 
solution for the LDF-A. To illustrate the procedure the LDF-A will be con¬ 
structed explicitly for example LDF’s of three simulated showers in Section 
10. The predictions of the constructed LDF-A will be compared with the 
polar density of the simulated showers. In Section 11 the paper is concluded 
with a brief summary. 

2 Cylinder model 

Atmospheric attenuation has a large effect on inclined showers |3]. One of 
the consequences is that it shifts the center of an elliptic iso-density contour. 
To model it we will hrst assume that all the particles run parallel with the 
shower core at the moment of arrival. Furthermore, we assume that contours 
of equal density are circles in the plane perpendicular to the shower direction. 
For the coordinate system we take the x and y axes in the horizontal plane 
and the z axis in the upward vertical direction. The origin is taken at the 
position where the shower core axis intersects the horizontal plane. The 
azimuth angle of the shower is, anti-clockwise, with respect to the positive 
x-axis. Without loss of generality we consider inclined showers with zero 
azimuthal angle, thus with the shower core in the x,^;-plane. The situation 
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Figure 3: Front of inclined air shower falling on a horizontal surface. 


is schematically shown in Fig. The tilted circle is perpendicular to the 
shower direction. At the moment the core reaches the surface in the origin 
of the coordinate system the shower front intersects the horizontal plane at 
the y-axis. We take a point N on the tilted circle. Its distance with respect 
to the origin is r. If the direction of the shower particles is conveniently 
assumed parallel to the shower core, the projection of N on the horizontal 
plane is point P. M is a point on the y-axis with identical y-coordinate as 
N and P. The angle between NM and PM and the angle between NP and 
the vertical axis both are equal to the zenith angle 9. From the geometry it 
follows 


MN‘^ = ON'^ - 0M‘^ = r‘^ -y^ J 

(1) 

2 2 2 /3,2 
r = X COS 9 + y 

(2) 

NP^ = x^ sin^ 9 

(3) 


Here and in the sequel (x,y) denote the coordinates of P in the horizontal 
plane. 

Without attenuation the asymmetry would be solely caused by the projection 
of the shower plane onto the horizontal observation plane. Alternatively, the 
intersection of a slant cylinder with a horizontal plane is an ellipse. The 
projection along the shower core axis means that the density along NM is 
projected to the larger PM. As a consequence the density v at the horizontal 
plane is smaller than the density p of the inclined shower front by a factor 
cos 9: 

v{x, y) = p{r) cos 9 . (4) 
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At the same time the iso-density contours are stretched to ellipses satisfying 
Eq. ([^. The horizontal ellipse and the inclined circle intersect each other 
and the positive y-axis at K. Denoting the y-coordinate oi K as k we obtain 

cos^ 9 + y^ = . (5) 

This is an ellipse whose semi-major axis a and semi-minor axis b are related 
to each other via b = a cos 9 and where the center of the ellipse coincides 
with the shower core. 

Denoting the x-coordinate of the shifted center as xm the general equation 
for a shifted ellipse is 


(x — xm)^ cos^ 9 + y^ = b^ , (6) 

where b is the size of the semi-minor axis. Since y = A; if x = 0 we also 
obtain 

x\j cos^ 9 + k^ = b^ . (7) 

From the latter two equations we can write the equation of the shifted ellipse 
also as 

(x^ — 2xxm) cos^ 9 + ip' = k‘^ . (8) 

The equation can be solved for k after we have determined xm as a function 
of k and 9. By means of the solution the LDF-A can be constructed. 

3 Modeling attenuation 

At the early stages of the longitudinal development the size of a shower 
increases. After the shower size has reached a maximum it approximately 
falls of exponentially with atmospheric depth. The attenuation length A is 
about 185 g cm’ A consequence of the attenuation of the shower 

during the traverse from N to P is that the density of shower particles 
is decreased by a factor where AX is the additional atmospheric 

depth met by shower particles between N and P. The atmospheric depth 
exponentially decreases with altitude with a characteristic length of about 8 
km. Except for shower with a very large energy and very large inclination 
the atmospheric depth between N and P approximately is constant. At the 
surface of the earth the increase AX is approximately equal to 0.13 g cm“^ 
for every meter travelled through the air. Hence, 

i/(x, y) = p{r) ■ ■ cos 9 , (9) 

where ^ = 0.13/185 7.0 • 10“'^ m“^. With the substitution of the Eq. ^ 

for NP this is: 

z/(x,y) =/9(r) • • cos0 , (10) 
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where, according to Eq. ([^, r = cos^ 9 + y'^. This is the basic equation 
for the analysis. It accounts for the attenuation at the late part of the inclined 
shower and for the reverse at the early part. As a consequence it leads to a 
shift of the elliptic density in the horizontal plane. The performance of the 
reconstruction of shower core positions should improve if the polar density 
function is modihed for the shift. In our coordinate system the late part of 
the shower is at the negative x-axis. Notice that a negative value for x leads 
to a decrease of the density. For x = 0, y = k we have 

z^(0, k) = p{k) ■ cos 6 . (11) 

To obtain the iso-density contour in the horizontal plane through (0, k) we 
have to solve the equation v{x,y) = u{0,k). That is, we have to solve the 
equation 

h'{x,y) = p{k) ■ cos 9 (12) 

or, more explicitly, 

p{r) ■ = p{k) , (13) 

where r = a/x^ cos^ 6 + y‘^. In the next section we will derive analytically a 
hrst order solution for this equation. 


4 Analytical approximation 


The key in the following analysis is the observation that the lateral density 
can be roughly described by the following exponential function: 

p{r) oc . (14) 


In Fig. I^the polar averaged combined density, binned with bin-width 1 m, 
and their approximation by the exponential function are plotted for three 
showers: one with energy 10^® and zenith angle 30°, one with energy 10^^ 
and zenith angle 45° and one with energy 10^® and zenith angle 52.5° re¬ 
spectively denoted as shower a, b and c. 

It is evident that the approximation is not particularly good. The expo¬ 
nential function is solely intended as a toy function. The values for the 
parameters of shower a, b and c respectively are 0.022, 0.026, and 0.017 
meter for ro and 0.25, 0.25, and 0.22 for w. In the remainder of the analysis 
we will solely use ro = 0.025 m and w = 0.25. By means of the ‘toy’ function 
for the lateral density the equation zz(x, y) = i^(0, k) can be written as 

g-(r/ro)“gC3;sine _ g-(fc/ro)™ ^ 


Hence, 


r^ = k^ 


/ Cxsin0\2/“' 

V (k/ro)^J 


(16) 
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Figure 4: Polar averaged lateral density of electrons and muons together for 
the three showers a, b and c as given in the text. The dashed curves are the 
approximations by the exponential function as given in the text. 


The latter can be expressed as 

{1 + xwT] cos , 

where 


h = 




tan 9 


(17) 


(18) 


w{k/rQY 

Since xwrj cos 9 < 1 for x smaller than 10^ m we will take a hrst order 
approximation for the right hand side of Eq. 0: 

PS(1 + 2xr/cos0) . (19) 

With the substitution of Eq. (i) for r and some rearrangement we obtain 


{x cos 9 —k^T]^ +'lY ^ k^ [l + k^rY^ . 
From the comparison with Eqs. (|^ and 0 we hnd 




and 


xm 


k'^rj 
cos 9 


( 20 ) 


( 21 ) 


( 22 ) 


With the substitution of the Eq. (18) for rj we obtain the following model 
prediction for the shift: 


^•^0 i2-w tan(9 
w cos 9 


(23) 
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Since xm > 0 the shower attenuation does shift the center of the ellipse 
towards the early part of the shower. The shower attenuation did not change 
the eccentricity, which is equal to sin0. Substituting ^ = 7.0 • 10“^ m“^, 
ro = 0.025 m and w = 0.25, we obtain as the ‘cylinder’ model prediction for 
the shift in m: 

o 1 'rr tan 0 . . 

XM = 1.1 • 10“^ • -^ . (24 

costl 


5 Cone model 

Another model is the cone model. That is, we consider paths from apex A 
through the horizontal plane at ground level. The differences in experienced 
atmospheric depth, due to the different path lengths, will be translated in 
attenuation. To this end we consider a shower cone with apex A in the 
same coordinate system as in Fig. As before, we regard inclined showers 
with zero azimuthal angle, thus with the shower core in the x,z-plane. For 
positions inside the cone the opening angle is denoted as 5. The situation is 
schematically shown in Fig. 

A 



Figure 5: Front of inclined air shower falling on a horizontal surface. 


As for the cylinder model, the tilted circle is perpendicular to the shower 
direction. We take a point N on the tilted circle as shown in Fig. Its 
distance with respect to the origin is denoted as r. The projection of N, 
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along the cone, on the horizontal plane is point P. M is a point on the 
y-axis with identical y-coordinate as N. The angle between NM and QM 
is equal to the zenith angle 9. The angle between NP and the vertical axis 
is equal to 0 + 5. From the geometry it follows 


{xA, VA, za) = {L sin 9, 0, L cos 9) 


(25) 


and 


{xn-, Un, zn) = T tan S (— sin/3 cos 9, cos (5, sin/3sin6l) , 


(26) 


where 


sin [3 


MN 

~ON 


and ON = r. The line through A and N intersects the z 
P with {x,y) coordinates 


(27) 


0 plane in point 


{x,y) 


Ltan5 

1 — tan 6 sin /3 tan 9 



(28) 


As before we let k be the y-coordinate where the projected contour intersects 
de positive y-axis. Using the coordinates as given before we find, to first order 
in 5, the following approximate lengths of paths AP and AK: 


AP ~ L(1 + (5 tan 9 sin /3) 


(29) 


and 

AK . (30) 

Taking 1030 g cm“^ for the atmospheric depth at ground level the slant 
atmospheric depth experienced in these paths is 


Nap 


1030 • (1 + (5 tan 9 sin /?) 


and 


Nak 


cos 9 

1030 
cos 9 


(31) 


(32) 


The difference AX between the atmospheric depth experienced by path AP 
and path AK is 


AX = Xap - Xax 


1030 • (5 tan 0 sin/3 
cos 9 


(33) 


This corresponds to an additional attenuation given by 

, (34) 

which can be elaborated to 

(35) 


g—AX/A _ g—7.9-10®^(5 sin/Stan 0 cos ^9 
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For S we have 


r r cos 9 


where h is the altitude of the apex in m. Hence 

g—AX/A g—sin/3 tan 6 

Since rsin^S = MN ~ —xcosO the attenuation can be written as 

g-AX/A _ g7.9-103/i-l5a;sin6» 

Equating v{x^y) with Ty{{),k) we obtain 

g-(r/ro)™g7.9'103/i“i^a:sin6» _ ^-{k/ro)^ 


(36) 


(37) 


(38) 


(39) 


The exponent of the attenuation differs only by a factor 7.9 • 10^from the 
one in the previous model. Proceeding in a similar way as in the previous 
section will therefore lead to a shift which is 7.9 • 10^times as large as 
the one of the previous model: 


7.9-10^ ^ , 2 j/i tan0 

XM = -7- --k^ -X 

n vw cos 9 


(40) 


For a shower with zenith angle 60° it is found for the difference in atmospheric 
depth between the late and early part of the shower 2AX = 370 g cm“^ at a 


distance of 1000 m from the core |7|- From equations (|33[) and (36) it follows 


for the difference between the late {j3 = 7^12) and early (/3 = —'xl2) part: 


2AX 


2 • 1030 • r tan0 
h 


(41) 


For 2AA = 370 g cm“^, r = 1000 m and 9 = 60° the latter equation is 
satished if h = 9.6 • 10^ m. If we substitute h = 9.6 • 10^ m, ro = 0.025 m“^, 
w = 0.25 and ^ = 7.0 • 10“^ m“^, the ‘cone’ model prediction for the shift in 


m IS 


xm = 9.2 


10-4 . ^1.75 . ^ 
COS 9 


(42) 


The latter is only 16% smaller in comparison to the cylinder model predic¬ 
tion. It can be imagined that the energetic particles near the core point on 
average to an apex with larger. That would correspond to an even smaller 
prediction. We return to this in Section 7 where we discuss the situation for 
h being a function of k. Anyway, in the following we will compare the shift 


solely with the cylinder model prediction (24), which is equal to the cone 
model prediction (40) if h ~ 7.9 • 10^ m. 
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6 Method 


In this section we describe the method of investigation of the shift in the 
lateral density. The method can be applied for the electron density, muon 
density and the combined density of electrons and muons together. We 
restricted ourselves to lateral densities of proton initiated showers. The 
showers were generated with CORSIKA-v7.4 [^, with QGSJET-II-04 j9l 110) 
+ GHEISHA CH for the hadronic interactions. The showers were generated 
without thinning. The horizontal observation level was set to 10 m. The 
energy cuts are 0.3 GeV for hadrons and muons and 3 MeV for electrons 
and photons. For each shower the lateral distribution is binned with bin size 
equal to 10/y^ m. As an illustration the binned lateral combined density of 
an arbitrary 10^^ eV shower with 45° zenith angle and 0° azimuth angle is 
shown in Fig. 



—400 —200 O 200 400 

X [m] 


Figure 6: The asymmetric lateral density of a 10^^ eV proton shower with 
zenith angle 45° and azimuth angle 0°. 


From the binned density, smoothened by means of a Gaussian filter with one 
bin (x X x) as sigma, the iso-density contours are determined. By means of 
minimization of the sum of squares the contours are fitted by an ellipse with 


12 


equation 


2 


(43) 



where a and b are the semi-major and semi-minor axes respectively. For the 
example shower of Fig. the hnal contour with density < p >= 1 m“^ 
is shown in Fig. [^together with the ellipse resulting from the ht. In this 
way we obtain values for the semi-major axis, the semi-minor axis and the 
value of xm- The center of the ellipse is denoted as M. The shower core 
is at the origin O. The focal points Fi and F 2 are at distance c from the 
center M. This distance is related to the semi-major and semi-minor axis 
via — 6^. 



-300 -200 -100 0 100 200 300 

X [m] 


Figure 7: The iso-density contour for density p = 1 m“^ of a 10^'^ eV proton 
shower with zenith angle 45° and 0° azimuth angle (solid). It is with good 
approximation equal to an ellipse (dashed). 


Next to Xm the ht also delivers a value for the pM coordinate. Its value, 
which is close to zero as it should, will be left out of the analysis. Ignoring 
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Um the Eq. (43) can be written either as 


or as 


(^{x-xm)^ +y^ = h^ 


(44) 



2x ■ xm 



(45) 


where 

= b^ - ^xii . (46) 

So, having determined the semi-major axis a, the semi-minor axis b and the 
shift Xm, we also know the corresponding value of k. 


From the model analyses, we expect b = a cos 9. In Fig. the value of 
6/a is plotted against zenith angle for the combined density of the simulated 
showers. 



^ (°) 


Figure 8: The ratio of the semi-minor and semi-major axis of elliptic lat¬ 
eral combined density of simulated showers for different zenith angles. The 
dashed line is cos 9. 


We see b/a follows cos 9 for small zenith angles. For larger zenith angles a 
small deviation shows up. The deviation increases with zenith angle. Since 
the expectation b = a cos 9 is based on polar symmetric iso-density contours 
in the front plane, the deviation suggests the iso-density contours in the front 
plane to be slightly elliptic with the major axis perpendicular to the azimuth 
direction. With the substitution of cos9 for b/a the Eqs. (44) through (46) 
reduce to the Eqs. (|^ through ([^. 


The method has been tested with artificial Poisson randomized shifted ellip¬ 
tic densities in order to check if the imposed shift is returned. The deviations 
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between the imposed and returned shifts were small, around 1 m for den¬ 
sity of 1 m“^ or less. Next to the inaccuracy of the method there also are 
contributions to the deviations due to the fluctuations of the densities. As 
a measure for the uncertainty the deviations of i/m with respect to the ex¬ 
pected value 0 are taken. For the simulated showers we found the standard 
deviation of i/m to depend on density roughly as ay ~ Assuming the 

standard deviation of xm to be comparable, we use it for the size of the error 
bars in the diagrams in the next section. 


7 General Monte Carlo results 


In this section we consider the lateral densities of two proton initiated show¬ 
ers with energy 100 PeV and 10 PeV both with zenith angle 45°. For both 
showers the shifts were determined for the electron density, the muon density 
and the combined density. For the combined density and the muon density 
the shift was determined for densities 0.001, 0.002, 0.003, 0.004, 0.006, 0.008, 
0.01, 0.02, 0.03, 0.04, 0.06, 0.08, 0.10, 0.20, 0.30, 0.40, 0.50, 0.64, 0.81, 1.00, 
1.44, 2.00 and 5.0 m“^. For the electron density the same densities were 
used with the densities 0.0002, 0.0003, 0.0004, 0.0006 and 0.0008 added 
to that. The shifts are shown in Fig. [^respectively Fig. 
curve in Figs, j^and 10 is the model prediction 


10 The dashed 


Comparing the model prediction with the determined shifts of the electron 
density we see the model prediction follows to a certain extent the shifts as 
determined for the electron density. The model predicts too low for k < 700 
and k < 500 m for the 100 PeV reps. 10 PeV shower. 


There are many reasons for the model to deviate from the determined shifts. 
To begin with, the plotted model prediction for the shift was based on a 
constant value for h. In reality h will depend on the distance to the shower 
core, as visualized in Fig. 2 of m- This suggests a large value for h near 
the core and a decreasing value for h for increasing k. According to the cone 
model a larger value for h implies a smaller value for the shift. As a con¬ 
sequence it will enhance the underestimation near the core. A small effect 
has the atmospheric depth decreasing exponentially with altitude. This will 
flatten the model curve for large k. We just took a constant value for the 
atmospheric depth in the model which is sufficient for our region of interest: 
k < 1000 m and 9 < 60°. 

The most important reason for the bad prediction probably is the local varia¬ 
tion of the attenuation. It can be imagined that the attenuation is large near 
the core and decreases for increasing k. The latter would enhance the shift 
near the core and flatten the model curve further away from the core. Al¬ 
ternatively, it might bring the model curve more in agreement with the shift 
curve of the combined density. To model it requires the knowledge of ^(fc) 
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k [m] 


Figure 9: The shifts of the electron, muon and combined density of a 100 
PeV shower with zenith angle 45°. The data points are connected with line 
pieces to guide the eye. The dashed curve is the model prediction. The 
densities 0.01, 0.001 and 0.0002 are depicted. 
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k [m] 


Figure 10: The shifts of the electron (triangle), muon (diamond) and com¬ 
bined (circle) density of a 10 PeV shower with zenith angle 45°. The dashed 
curve is the model prediction. The numbers 0.01, 0.001 and 0.0002 depict 
the density. 
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as a function of k. The local h{k) can possibly be obtained from simulated 
showers by inspection of the directions of the electrons when they arrive at 
the observation plane. The knowledge of ^(fc) seems more difficult: besides 
the inspection of the local energy distribution of electrons it also requires a 
relation between the distributions and the local attenuation. On the other 
hand, if one succeeds in describing and modeling the atmospheric depth and 
h as functions of k, parameterized by zenith angle, there is an opportunity to 
retrieve the local attenuation ^(k) from the shift XMik) as determined from 
the simulated electron density. For a model which predicts the shifts of the 
muon density one has to consider the decay of muons to electrons and the 
subsequent atmospheric attenuation of the electrons. The combined density 
the shift then follows from 


XM,e+^ {k') 


PeXM,e{^) + PtJ-XM,pi{k) 

Pe{k)+ p^i{k) 


(47) 


Because of the aforementioned reasons it is difficult to derive a precise model 
for the shifts of the electron density, let alone for the muon density and the 
combined density. Therefore we will not proceed in that direction. Instead, 
we will focus our attention on the behavior of the determined shift for rela¬ 
tively large densities. 


In Fig. 11 the shift curves of the combined density of a 10 PeV and a 100 
PeV shower with zenith angle 45° are once more plotted. They are identical 
to the ones in Fig. and Fig. except that the dots and error bars are 
left. 

The slope of both curves show some curious irregularities. For the 100 PeV 
curve these are around k = 200, 580 and 1400 m. For the 10 PeV curve 
we see them around k = 250 and 650 m. The question arises whether these 
irregularities are the remnants of consecutive hadronic interactions. 


The two shift curves fall on top of each other for k < 200 m. At k = 200 m 
they do branch. Beyond the fork the difference between the curves slightly 
increase for increasing k. In general this means that the shift curves are 
not independent of shower size. The density of the 10 PeV shower at the 
branching point is about 0.06 m“^. In Fig. 
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the shifts of the combined 
density as found for the 10 PeV shower is plotted on top of the ones for the 
100 PeV shower. For the 10 PeV shower the plotted densities (white) are 
0.2, 0.3, 0.4, 0.5, 0.64, 0.81, 1.00, 1.44, 2.0, 5.0 and 10 m-^. For the 100 PeV 
shower the plotted densities (black) are 0.3, 0.4, 0.5, 0.64, 0.81, 1.00, 1.44, 
2.0, 5.0, 10, 20 and 50 m-^. 

We see the shift curves practically fall on top of each other within the given 
density domains. This means that we can try to find a relation between xm 
and k independent of shower size (or energy) similar to the model prediction. 
In addition, the curves are almost linear. As we will see further on the latter 
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allows for an analytical solution for the shifted polar density. 

8 Specific Monte Carlo results 

In this section we will investigate the shift in combined densities of a set of 
simulated showers. The energies of the showers are 10^^, 10^®, 10^^ and 10^® 
eV. The zenith angles of the showers range from 7.5° through 60°, in steps 
of 7.5°. As an illustration the ratio Ne/(Ne + is plotted against energy 
for several zenith angles in Fig. The energy - zenith angle entries are 
shown as black dots. 



Logio(E[eV]) 


Figure 13: The ratio Ne/{Ne + against energy for zenith angles 0° 
through 60°, in steps of 7.5°. 


Not all possible combinations of energies and zenith angles are applicable 
for investigation. For showers with relatively low energy and relatively large 
zenith angle, shown as the gray region in Fig. 13 the small shower size at 


observation level does in general not allow for a determination of the shift. 
The choice for generating showers without thinning is made to avoid possi¬ 
ble deviations caused by thinning. The consumption of computer time and 
of storage space grows exponentially with the size of the simulated shower 
|13) . This is extremely the case for shower simulation without thinning. As 
a consequence the library of showers generated without thinning is limited, 
in particular for large energies. For the largest energy considered, 10^® eV, 
the library is momentarily limited to 10 showers for zenith angle 60°, 10 
for zenith angle 52.5°, 8 for zenith angle 45°, 5 for zenith angle 37.5° and 
none for zenith angles 30°, 22.5° and 15°. To obtain a sort of equal share 
in our diagrams we take 10 showers for each of the other energy - zenith 
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angle entries. For each shower we determine the iso-density contours for 
combined densities 0.50, 0.64, 0.81, 1.00, 1.44, 2.0, 5.0, 10, 20 and 50 m“^ 
for as far as these densities occur in a shower, thus maximum 10 data points 
per shower. For the densities considered this is close to the shift as we would 
have obtained it from the electron density, except for regions were the muon 
component dominates: for large distances to the core and for zenith angles 
in the neighborhood of 60° and larger. As depicted in Fig. 13 for zenith an¬ 
gle 7.5° through 30° showers were used with energy 10^^, 10^® and 10^^, for 
zenith angle 37.5° through 52.5° showers were used with energy 10^®, 10^'^ 
and 10^® eV and for zenith angle 60° showers were used with energy 10^'^ 
and 10^® eV. For 10 showers at 3 energy decades we obtained a maximum of 
300 data points for each zenith angle. 


As we will see, and as suggested by the model results, for each zenith angle 
the data points follow a curve independent of energy. The curves are htted 
with a function similar to the one resulting from the models. To be specihc, 
for each zenith angle we plot the xm against k and ht the result by the 
equation 


Xm = 


B 


taxi 9 
cos 9 


(48) 


As an illustration the Xm are plotted against k for zenith angle 30° and htted 


with Eq. (48), see Fig. 14 



k [m] 


Figure 14: The shift xm against k (unhlled circles) and the best ht curve for 
zenith angle 30°. 
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In the latter figure the error bars since they are in most of the cases smaller 
than the size of the plot markers. The xm grows with k along a curve in¬ 
dependent of the primary energy of the showers. For the 1 PeV showers the 
data points are in the region fc < 80 m. For the 10 PeV and 100 PeV the 
regions are k < 250 m and k < 500 m respectively. For the parameters of 
the ht we hnd A = 0.020 and B = 1.28. For the goodness of ht we hnd for 
the Pearson test 0.98 as the p-value. These hgures hold for zenith angle 
30°. For zenith angle 7.5° through 30° the diagrams are shown in Fig. 

For zenith angle 37.5° through 60° the diagrams are shown in Fig. 16 



k [m] 




k [m] 



k [m] 


k [m] 


Figure 15: The shift xm against k (unhlled circles), the best ht curve (solid), 
the model prediction (dotted) and a linear curve (dashed) for zenith angle 
7.5° through 30°. 


In both diagrams the model predictions for the shift of the electron densities 
are plotted as well for reasons of comparison; for small zenith angles the 
combined distribution is dominated by electrons. 

For small zenith angles, 7.5° and 15°, the spread of the shifts are mainly 
governed by the uncertainty of the measurement. For large zenith angles. 
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80 



Figure 16: The shift xm against k (unfilled circles), the best fit curve (solid) 
and the model prediction (dotted) and a linear curve (dashed) for zenith 
angle 37.5° through 60°. 


e 

A 

B 


p-value 

datasize 

7.5° 

0.010 

1.41 

25.7 

0.11 

248 

15° 

0.013 

1.36 

8.74 

0.95 

220 

22.5° 

0.019 

1.29 

11.7 

0.86 

274 

O 

O 

CO 

0.020 

1.28 

7.79 

0.98 

253 

37.5° 

0.037 

1.17 

12.5 

0.82 

250 

45° 

0.043 

1.13 

10.0 

0.93 

267 

52.5° 

0.034 

1.13 

9.21 

0.91 

199 

O 

O 

0.019 

1.09 

11.1 

0.74 

156 


Table 1: Parameters A and B of the fits, the and p-values and the number 

of inspected contours for zenith angles as given in the first column. 
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52.5° and 60°, the spread of the shifts are mainly due to shower to shower 
variations. For each zenith angle the values A, B, the p-value have 

been tabulated, see Table 


The values of B as given in Table decrease for increasing zenith angle. 
The values of B being close to unity suggests to consider a linear relation 
between Xm and k. From hts with 

XM = C ■ k - (49) 

cosU 


it is found that C scales as cos0. On average C ~ 0.116 cos 0, except for 
6 = 60° where C ~ 0.075 cos 0. Writing C as 25 cos 0, the proposed linear 
relation is as follows 


Xm = 25 • /ctan0 , 


(50) 


where 5 = 0.058. For 0 = 60° 
linear relation is shown in Figs. 


the value of 5 is about 35% smaller. The 
15 and 16 as a dashed curve. For small 
zenith angles the linear approximation overestimates the shift for A; ~ 150 
m, the difference being just about 1 m. For zenith angle 37.5° and 45° 
it underestimates by 20% in the region where the density is small, p < 1 
m“^. Taking the inaccuracies for granted, a linear equation is advantageous 
since, as we will see further on, it allows for an analytical solution for the 
description of the LDF-A. At the end of Section 10 a remark will be made 
about the possible application of the more accurate power law (48). 


9 The polar density 


In this section we will perform the conversion of an LDF to an LDF-A. 
Substitution of the shift (50) in Eq. (l8|) gives 


(^x^ — Ax ■ S ■ k tan 0^ cos^ 9 + y^ = k^ . 
Solving for k we obtain 


(51) 


k = —Sx sin(20) + + x‘^{l + 45^ sin^ 0) cos^ 0 , (52) 

For 5 = 0.058 the term 45^ sin^ 0 is negligible with respect to 1. With good 
approximation we therefore have 


k = —5xsin(20) -|- \J+ x'^ cos^ 0 


(53) 


According to Eq. (12) the polar density z^(x, y) including the shift is obtained 
in Cartesian coordinates by substituting expression (53) for k in the polar 
symmetric density p and by multiplying it by cos0. The horizontal polar 
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density v can also be written in polar coordinates. With the substitution of 
X = rcosa and y = rsina the Eqs. (51) and (53) respectively read 

cos^ a cos^ 9 — 2r ■ S ■ k cos a sin(20) + sin^ a = (54) 

and 

k = —Sr cos a sin(20) + r\Jl — cos^ a sin^ 6 . (55) 

The second term, the square root part, is due to the ellipticity of the density 
as caused by the projection M- The hrst term on the right hand side of Eq. 
(55) is due to the shift. The polar density including the shift is obtained in 
horizontal polar coordinates in the same way as for Cartesian coordinates: 

zz(r, a) = p{k) cos 9 (56) 


with k as given by (55). 


To obtain the polar density v we need the polar symmetric density p. A 
good approximation for p is found by polar averaging the horizontal den¬ 
sity and htting it with a suitable LDF. For electromagnetic showers a well 
known LDF is the one of Nishimura, Kamata and Greisen (NKG) |15[ I16| . 
Most LDF’s are modihcations of the NKG function HZl I18| . For muons a 
well known lateral density function is the one of Vernov |19| . However, it 
can also be described by a NKG type of function 1201 [2T| . For radii smaller 
than about 300 m the combined density of electrons and muons can also be 
described by an NKG type of LDF: 


Pe+ix{r) = Ne+f, • c - /(r) , 


(57) 


where 



(58) 


is the structure function and where c usually is the normalization. For¬ 
mally this LDF is similar to the one used for the KASCADE experiment 
|18) . There the quantities s — a and s — f3, with s the shape parameter 
(a remnant of the age parameter), play a similar role as si respectively S 2 - 
The parameter ro plays a similar role as the Moliere radius in the original 
NKG function. From the simulated showers it is found that tq is close to 30 
m. Fixing tq to 30 has only a marginally effect on the ht values for si and S 2 - 


For radii larger than about 300 m it underestimates the combined density. 
The deviation is caused by the relatively large muon component. To adjust 
for the muon component we let us motivate by the Greisen function j22| . 
That is, we multiply the LDF by (1 -|- xyjy^)- The latter multiplication 
complicates the normalization. We therefore take ro = 30 and use c as a ht 
parameter like si and S 2 - Thus 
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10 Comparison with simulated densities 


The performance will be illustrated by means of the same three showers a, b 
and c as already used in Section 4. Their polar averaged combined densities 
and their fitting curves are plotted in Fig. As before the polar averaged 
densities are binned with bin-width 1 m. 



Figure 17: Polar averaged combined densities for the three showers a, b and 
c as given in the text. The dashed curves are the fits with the LDF as given 
in the text. 


We see the fit curves follow the combined densities also beyond 300 m. The 
number of electrons and muons and the values found for the parameters are 
shown in Tabled 


shower 

e 

^e+ii 

c 

Si 

S2 

a 

CO 

o 

o 

490538 

0.000283 

-0.441 

-2.787 

b 

45° 

1967956 

0.000205 

-0.355 

-2.638 

c 

52.5° 

5834026 

0.0000903 

-0.309 

-2.253 


Table 2: The zenith angle, the shower size and the fit values for c, si and S 2 
for the combined density of the showers a, b and c as given in the text. 


To obtain the LDF-A of the three showers we multiply the polar averaged 
LDF with COS0 and replace r for k. The final prediction for the LDF-A 
including the shift is 

= JV.+. ■ “'O ■ ■ (s) ' + s) ’ 
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where k is as given by (55) and with the parameter values as given in Table 


For the three simulated showers we inspected the polar variation of the com¬ 
bined density at different radii. To this end the density was binned with bin 
size 7rr/12 in the angular direction. For the bin size in the radial direction 
we took 1, 1, 2, 3, 4 and 5 m for radii 10, 20, 50, 100, 200 and 500 m re¬ 
spectively. The bins were just large enough to balance out to some extent 
the Poisson fluctuations in the density. For shower a, b and c the result is 
plotted in Fig. [T^ through |20| together with the LDF-A prediction. 



Figure 18: The binned combined density of shower a (10^® eV , zenith angle 
30°) against polar angle for different radii (dots) and the constructed LDF-A 
prediction (solid). 


We see the constructed LDF-A nicely follows the polar density. Angular in¬ 


dependent deviations, such as in the upper left panel of Fig. 19, are caused 
by the inaccuracy of p. Any inaccuracy in the underlying LDF will be re¬ 
flected in the LDF-A. This does not take away that an LDF-A including the 
shift still follows the polar variation of the density better than an LDF-A 
without the shift. 


The accuracy of the LDF-A for the three example showers suggests that 
it is probably sufficient to consider a single linear relation between xm and k 
independent of zenith angle. If a better accuracy is desired one has to apply 
the power law (48) with ht values for A and B as given in Table The latter 
approach requires the ht coefficients either to be tabulated for different 6 or 
to be parameterized to 6 by means of a suitable function. Next to this, for a 
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-180 -90 0 90 180 


a(°) 

Figure 19: The binned combined density of shower b (10^^ eV , zenith angle 
45°) against polar angle for different radii (dots) and the constructed LDF-A 
prediction (solid). 





Figure 20: The binned combined density of shower c (10^® eV , zenith angle 
52.5°) against polar angle for different radii (dots) and the constructed LDF- 
A prediction (solid). 
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shift as given by the power iaw the equation 54 wouid have been as foiiows 
cos^ a cos^ 9 — 2r cos a ■ A - sin 0 + sin^ a = . (61) 

For B a non-integer the iatter equation has to be soived numericaiiy for k. 


11 Summary 


For eiectron densities the shift of the center of eiiiptic iso-density contours is 
modeied. For combined densities of simuiated showers the shift is determined 
for different zenith angies. An approximate iinear reiation between the shift 
and k aiiows for an anaiyticai soiution for the asymmetric iaterai density. 
The conversion of an LDF to an LDF-A inciuding the shift is described. The 
conversion consists in two steps corresponding to the effects of the projection 
and of the attenuation. The hrst step is the multiplication of the LDF by 
cos 0, where 6 is the zenith angle. The second step is to replace r by k, where 
k is given by Eq. (55). The result is an LDF-A for the situation where the 
azimuth angle is equal to zero. The LDF-A for a non-zero azimuth angle (p 
requires a third, trivial step: replacing the polar angle a hy a — p. 


Zenith angle 60° is a sort of transition point. Below this point there is a 
relatively large shift while the effect of the geomagnetic held is negligible. 
Above the transition point the muon component becomes dominant. As a 
consequence the relative shift is smaller. At the same time the inhuence of 
the geomagnetic held rapidly grows with zenith angle. Above the transition 
point the consequences of the shift for the asymmetry of the density will 
be overwhelmed by the effect of the geomagnetic held. The inhuence of the 
geomagnetic held on very inclined air showers requires a different modeling 

[71IM1IM1ES]. 


The aim of the paper was to consider the situation for zenith angles smaller 
than 60°, where the shift is mainly governed by the attenuation of the elec¬ 
tron component and therefore to a certain extend substantial. The inclusion 
of the shift leads to a more accurate description of an asymmetric polar 
density. It therefore may be worthwhile to take the shift into account for 
reconstruction purposes. It should be emphasized that the Eq. (50) and 
the shift part of Eq. (55) were derived for observation at sea level. For 
an observation level at a different altitude, the shift will be different. Even 
if the inclusion of the shift improves the accuracy of reconstruction only 
marginally, the contents of the paper may still contribute to the descrip¬ 
tion and understanding of horizontal polar densities of inclined cosmic air 
showers. 
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